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in the double-exchange spin ice model on the kagome lattice 
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We investigate a detail of a dodecamer cluster ordering in a double-exchange spin ice 
model on a kagome lattice. In frustrated systems, ordinary spin orderings are suppressed 
and macroscopic degeneracy remains down to low temperatures. In some frustrated systems, 
the degeneracy is lifted due to residual interactions and cluster orderings are stabilized. In 
the present model, the spin ice state is first formed at intermediate temperatures, and further 
entropies are released at lower temperatures as the dodecamer phase emerges. Since the spin 
symmetry is not broken in the dodecamer phase, there still exists macroscopic degeneracy. 
At further low temperatures, a possible spin ordering due to inter-dodecamer interactions is 
proposed. We discuss that such a multiple-site clustering larger than a bond-pair might be 
generic to frustrated systems where macroscopic degeneracy is lifted by residual interactions. 

KEYWORDS: dodecamer, cluster order, double-exchange model, kagome lattice, frustration, 
metal, spin ice, Monte Carlo simulation 



1. Introduction 

Intensive studies have been devoted to frustrated systems because of their unexpected 
behaviors. In spin systems, a frustration arises when all interaction energies between localized 
spins can not be minimized simultaneously in the classical picture. It is known that the 
frustration suppresses the long-range order and often induces the novel ground state. The 
macroscopically degenerate ground state accompanied with a residual entropy is one of the 
well-known states in frustrated Ising spin systems. Such a state is seen in the antiferromagnetic 
(AF) Ising model on a triangular, 1 a kagome 2,3 and a pyrochlore 4 lattices. On the other hand, 
a peculiar cluster ordered state can be induced by the frustration. In the Majundar-Ghosh 
model 5 and the Shastry-Sutherland model 6 which are frustrated quantum spin systems, the 
dimer-singlet states with a spin gap are realized as the unique ground states. The entropies 
in the systems are released by the dimer-singlet formations for all localized spins, which 
can be viewed as a cluster ordering since the systems are tiled by the dimer-singlets. Such 
dimer-singlet ground states are driven by geometrical frustrations and quantum effects of the 
interactions. The cluster ordering can be the way to stabilize the states in the frustrated 
systems. 

*E-mail address: shimomura@phys.aoyama.ac.jp 
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Fig. 1. (a) Dodccamer phase. Open and closed circles represent out- and in-spins for the up-trianglcs, 
respectively. The dodecamer is represented by the shadowed spin cluster. Each dodecamer is 
labeled i, j, k and others, (b) Twelve sites which construct the i-th dodecamer. 

Actually, such a cluster order is seen even in realistic materials. For example, the 
Shastry-Sutherland model is realized in the two-dimensional spin- 1/2 orthogonal dimer sys- 
tem SrCu2(BOs)2, where Cu sites are tiled by the dimer-singlet. This material is novel in 
terms of the existence of the exact dimer-singlet ground state with the finite spin gap in the 
two-dimensional system. 7 ' 8 Furthermore, as possible other experimental examples of the spin 
cluster order in frustrated spin systems, hexamers of Cr 3+ in a cubic spinel ZnCr2 04 9 and 
octamers of Ir 3+ and Ir 4+ in a thiospinel CuIr2S4 10 have been reported. 

In this way, cluster orders in frustrated spin systems are seen in both theoretical and 
experimental researches, and have attracted much attentions. Even in electron systems, similar 
cluster ordering induced by the frustration can be realized. As an example of such cluster 
orders in electron systems, we have reported dodecagonal localized spin cluster, "dodecamer" , 
order in a double-exchange spin ice (DESI) model on the kagome lattice 11 [see Fig. 1(a)]. In 
the DESI model, localized spins have local uniaxial anisotropies for triangles which constitute 
the kagome lattice: Localized spins are forced to point in the directions connecting the center 
of mass and the three vertexes of the triangles, that is, they point either inward (in-spin) or 
outward (out-spin) for up-triangles. There is an effective ferromagnetic interaction between 
localized spins due to the double-exchange (DE) mechanism. 12 The situation that there are the 
ferromagnetic interactions and the uniaxial anisotropies for localized spins is similar to that in 
the spin ice systems. 13 From an analogy with the spin ice systems, it is obvious that the DESI 
model has the frustration and a spin-ice-like behavior. The dodecamer formation is driven 
by the frustration and the kinetic energy gain due to the DE mechanism: At sufficiently 
low temperatures, electrons move along a certain selected path to gain the kinetic energy 
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under the frustration. The dodecamer formation is a consequence of the path selection by 
electrons in order to gain the kinetic energy. In the dodecamer phase, the kagome lattice can 
be completely tiled by dodecamers, which leads to the dodecamer order with the transnational 
symmetry breaking. Note that there are two kinds of dodecamers: One is the dodecamer with 
in-spins on the tips %2n-i = 1, 2, - - • ,6) shown in Fig. 1(b) and the other is out-spins. From 
the existence of two dodecamers, it is found that the spin symmetry is not broken in the 
dodecamer phase. 

In this paper, we focus on the behaviors of the entropy release in the DESI system. As 
the temperature is lowered, the thermodynamics in the DESI system is as follows. (1) The 
spin-ice-like state is stabilized at intermediate temperatures. (2) The dodecamer order appears 
at sufficiently low temperatures. In the spin-ice-like state, there still remain a large part of 
degrees of freedom for localized spins, i.e., there is a finite entropy which is almost consistent 
with the residual entropy of the AF Ising model on the kagome lattice. 3 A large part of the 
entropy is released by the dodecamer ordering. There is a possibility that the entropy still 
remains even in the dodecamer phase, since each dodecamer is two-fold degenerate. 

2. Double-Exchange Spin Ice Model 

In this section, let us introduce the DESI model on the kagome lattice, which is the 
Anderson-Hasegawa model 14 with local uniaxial anisotropics. The mechanism for the frustra- 
tion in this system is similar to that in the spin ice systems. 13 It is expected that similar 
degeneracy as in the spin ice systems exists in the DESI model. 
2.1 Anderson-Hasegawa model 

The Anderson-Hasegawa model 14 is a modified DE model with infinitely large Hund's-rule 
couplings: 



where c\{ci) is an electron operator which creates (annihilates) an electron with a spin parallel 
to a localized spin at site i, Si, and [i is the chemical potential in the grand canonical ensemble. 
Here, t(Si, Sj) is given in the form, 



where t is a transfer integral and 9y is a relative angle between the localized spins Si and Sj. 
The angle 9ij is given by 



where (9i,(pi) is an angle of Si in the polar coordinate. Note that t(Si,Sj) has a maximum 
value t when Si and Sj are parallel. Thus, ferromagnetic alignment of the localized spins 
broadens the bandwidth of the system, which produces a kinetic energy gain. Therefore, there 




(1) 




(2) 



cos Oij = cos &i cos 9j + sin 9i sin 9j cos(^j — 4>j)i 



(3) 
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(c) two relative angles (d) possible eight spin configurations 

Fig. 2. (a) Kagome lattice. The i-th unit cell surrounded by the dashed line consists of three sites 
i\, %2 and 13. a is a lattice constant, (b) Uniaxial anisotropics for localized spins. Localized spins 
have two degrees of freedom, "in" and "out", for a triangle, (c) two relative angles in the DESI 
system: The relative angle 9ij between in-spin and out-spin, and between the same spins, are ^ 
and respectively, (d) Possible eight configurations for an up-triangle in the kagome lattice. 



is an effective ferromagnetic interaction between the nearest-neighbor (n.n.) localized spins. 

We set the transfer integral t = 1, hereafter. 

2.2 Frustrated electron system : double- exchange spin ice model 

We consider the Anderson-Hasegawa model on the kagome lattice [see Fig. 2(a)], and 
introduce uniaxial anisotropies for the localized spins as follows: They are forced to be either 
in-spin or out-spin for up-triangles on the kagome lattice [see Fig. 2(b)]. Under the conditions, 
relative angles between the localized spins, 9ij (Eq.(2)), can take only two values, i.e., Q\a = 7r/3 
(the angle between an in-spin and an out-spin) and Of- = 2ir/3 (the angle between the same 
spins) [see Fig. 2(c)]: The transfer integral of electrons t(Si,Sj) has t\ = cos(6^-/2) = v3/2 
and t2 = cos(9?j/2) = 1/2. Thus, electrons can easily move through the ferromagnetic bonds 
with the transfer integral t\ due to the kinetic energy gain. 

To clarify the frustration in the DESI system, we introduce a pseudospin representation. 
We dehne a pseudospin corresponding to Si, 

+1 (S-i is an out-spin), 
Ti = { (4) 
-1 (Si is an in-spin). 
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In this pseudospin picture, the Hamiltonian (1) is rewritten in the form 



H = ^2 1 (n, Tj ) (c\cj + /i.e.) - fi^2 4°i> 

(iJ) 1 



(5) 



where 




(6) 



h (n = Tj). 



The effective ferromagnetic interaction between the n.n. localized spins is regarded as the 
effective AF interaction between the n.n. pseudospins. Thus, there is the n.n. AF interaction 
on the kagome lattice. Using the pseudospin picture, it is easily understood that the system 
has the frustration. 

2.3 Spin ice state in the double- exchange spin ice model 

In the DESI system, each triangle has four types of spin configurations, i.e., "three-in", 
"three-out", "two-in one-out" and "one-in two-out" as shown in Fig. 2(d). These correspond 
to "three-down" , "three-up" , "two-down one-up" and "one-down two-up" , respectively, in the 
pseudospin picture. Naturally, as the temperature is lowered, three-in and three-out configu- 
rations are unfavored by the ferromagnetic interaction due to the DE mechanism. In that case, 
the majority of localized spin configurations on the triangles are two-in one-out and one-in 
two-out. Here, each configuration is three- fold degenerate. A total spin state in the system is 
expected to have the macroscopic degeneracy since there still exist many degrees of freedom 
that each triangle satisfies either two-in one-out or one-in two-out. From the analogy with the 
spin ice systems, 13 an "ice state" in the DESI system is defined as a state with the macroscopic 
degeneracy and the two types of spin configurations, i.e., two-in one-out and one-in two-out. 

We can characterize the ice state by using a vector spin chirality V{ defined by 



where i\, i% and 13 represent three sites of the triangles in the kagome lattice shown in Fig. 
2(a). Here, we set \Si a \ = 1 (a = 1, 2 and 3). Vi has only z-component since localized spins in 
the DESI system are coplanar, where we consider that the kagome lattice is spread in xy-plane. 
vf has -3 (three-in and three-out) or 1 (one-in two-out and two-in one-out). It is expected that 
the thermal average of the vector spin chirality (v f) vanishes if the four types of configurations 
are realized equally. On the other hand, (vf) approaches unity when a large number of spin 
configurations for triangles are two-in one-out and one-in two-out. Furthermore, we define a 
fluctuation of vf, which is given in the form 



To see the existence of the ice state, we calculated the total vector chirality v z = ^2i(vf) /N^ 
and the total fluctuation Xz = Y2i x|/At- Here At represents the number of triangles on the 







(8) 



T 
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(a) T (b) T 

Fig. 3. (a) Temperature dependence of the specific heat, Cg, for the system with 4 x 4 (= 48 sites) 
(A) and 6 x 6 (= f 08 sites) (□) unit cells at fj, = 0. The specific heat has a double-peak structure, 
a broad peak (T ~ 0.2) and a sharp peak (T ~ 0.015). (b) Temperature dependence of the specific 
heat, Cg, and the entropy, Ss, (B) for the system with 6x6 unit cells. 

kagome lattice which is given by §As, where N$ is the number of sites. The summation runs 
over all up- and down-triangles on the kagome lattice. 

One may think that the DESI model on the kagome lattice could be completely mapped 
to the n.n. AF Ising model on the same lattice, which has a macroscopically degenerate 
ground state. 2 ' 3 However, the interaction of the DESI system are determined by the kinetics 
of electrons, and it is expected that the behavior of the DESI model is different from that of 
the frustrated Ising spin system with the n.n. interaction. 

3. Results : Thermodynamics of the Double-Exchange Spin Ice System on the 
Kagome Lattice 

In this section, we present the Monte Carlo (MC) results of the DESI model on the kagome 
lattice at finite temperatures. We applied a MC technique to the calculation of thermodynamic 
quantities by using the grand canonical ensemble [see Appendix]. We apply the Metropolis 
algorithm for the updates of spin configurations and typically run 100,000 MC steps for 
measurement after 10,000 thermalization steps. 

In ref. 11, we have already reported the dodecamer order of localized spins at sufficiently 
low temperatures. The dodecamer order is realized in a wide range of //, which is roughly 
estimated as —0.3 < fj, < 0.2. Thus, the DESI models around \i ~ are expected to show 
qualitatively similar behaviors. Therefore, we treat the DESI model at \i = 0, hereafter. 

Let us investigate the behaviors on the DESI model at finite temperatures. In order to see 
how the entropy of localized spins is released in the system, we have investigated a specific 
heat and localized spin configurations. As explained in Appendix, the specific heat of the 
system can be decomposed into an electronic part Ce(T) and a localized spin part Cs(T). 
The information of localized spin configurations is included in Cs(T). In particular, we are 
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interested in the behaviors of the localized spins at low temperatures. Thus, we only focus 
on Cs(T), which is defined as the specific heat per site hereafter. In Fig. 3(a), we show the 
temperature dependences of Cs(T), for the system with 4x4 unit cells (iVs = 48) and 6x6 
unit cells (N$ = 108) at /x = 0. Cs(T) has a double-peak structure: A broad peak (T ~ 0.2) 
and a sharp one (T ~ 0.015). In Fig. 3(b), the temperature dependence of the Cs and the 
entropy Ss obtained from Cs for the system with 6x6 unit cells are shown. These results 
provide us a picture that there exist two steps to release the entropy in the DESI system. 
3.1 Crossover to the ice state 

Let us consider the broad peak in Cg(T) at T ~ 0.2. As mentioned in Sec. 2.3, it is 
expected that the system exhibits an ice state as temperature is lowered. We speculate that 
the broad peak in Cg(T) is associated with the formation of the ice state. In order to clarify 
this idea, we investigate the entropy per site Ss(T) from Cs(T), 



Here we use S's(oo) = log 2, since all spin configurations with the number of states 2 s can be 
realized in the high temperature limit (T = oo). The data for S$(T) are shown in Fig. 3(b). 

In the high temperature region, Ss(T) deviates from log 2 at T < 1, since the interactions 
are mediated by the DE mechanism. In the intermediate temperature region, we see that 
Ss(T) has a plateau structure at T ~ 0.05 with Ss(T) ~ 0.5. Since the spin ice system on a 
kagome lattice can be mapped to an AF Ising model on the identical lattice, the entropy of 
the ice state for the DESI system is roughly estimated from the residual entropy of the AF 
Ising model on the kagome lattice, 3 Sising — 0.5, which is consistent with the observed value 
for Ss- 

An alternative clarification is given by the temperature dependence of v z and Xz- The 
results for the system with 6x6 unit cells (108 sites) and \x = are shown in Fig. 4. Xz has 
the maximum value at T ~ 0.2, and then, below which it rapidly decreases until T ~ 0.05. 
Similarly, it is found that v z has a inflection point at around T ~ 0.2, and saturates below 



These results suggest that the broad peak in Cs(T) at T ~ 0.2 is due to a crossover to the 
ice state. At T ~ 0.05, the ice state is stabilized, where a macroscopic number of thermally 
degenerate states still remain. From the analogy with the spin ice system, 13 it is considered 
that the origin of the ice state in the DESI system is the short-range (n.n.) interaction between 
localized spins. The ice state in which a large number of the entropy remains is indicated by 
the broad peak in the specific heat and is driven by the short-range ferromagnetic interaction. 
3.2 Transition to the dodecamer ordered state 

In this section, we focus on the sharp peak in the specific heat. The sharp peak means 
further entropy release in the system which results from the reduction of degrees of freedom 




(9) 



T ~ 



0.05. 



7/18 



J. Phys. Soc. Jpn. 



Full Paper 



1.0 
0.8 
^ 0.6 
0.4 
0.2 



□ 

x — Q 



~Z- 



6 
5 
4 

3 ^ 
2 
1 




0.0 0.1 0.2 0.3 0.4 0.5 0.6 

T 



Fig. 4. Temperature dependence of the vector spin chirality v z (d) and the susceptibility of the chi- 
rality XzC) f° r the system with 6x6 unit cells at \x = 0. i> 2 approaches unity as the temperature 
is lowered. \z has a broad peak at T ~ 0.25. 



for localized spins. As an origin of such an entropy release, it is considered that some kind 
of ordering of spins are generated due to the residual effects in addition to the short-range 
interaction in the model. An effective long-range interaction due to the kinetic energy gain 
of electrons is one candidate for such effects, since higher order hopping process of electrons 
is considered to be important at low temperatures. Indeed, the dodecamer order driven by 
the kinetic energy gain of electrons appears at low temperatures in the system. 11 Thus, it is 
expected that the sharp peak corresponds to a change from the ice state to the dodecamer 
ordered state. 

The order parameter for the dodecamer formation is characterized by using bond correla- 
tions, and is defined by 

12 6 \ 

^ ] T in T in+l + y 1 T «2m r i2m+2 J 1 (10) 

n=l m=l / 

where a site i n is indicated in Fig. 1(b) and we set 113 = i±, in = ii. Note that the order 
parameter di is a unity when the dodecamer is formed. The dodecamer structure factor is 
given in a form 

D q = ^J2 eiq<n ~ rj) ( d ^ (11) 

where Ac is the number of unit cells (Ac = ^Ag) and Yli j represents the summation for any 
pairs of the order parameter shown in Fig. 1(a). D q has three independent maximum peaks 
at Qi = *(0, 27r/v3a), Q2 = *(vr/o, 7r/\/3a) and Q3 = *(vr/a, — ir/y^a), those characterize 
the dodecamer order. An average of maximum peaks, Dq = {Dq 1 + Dq 2 + Dq 3 }/3, diverges 
increasing the system size at sufficiently low temperatures. 11 This result indicates that the 
dodecamer order is realized in the DESI system on the kagome lattice. Note that the change 
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Fig. 5. Temperature dependence of Dq and specific heat Cs for the system with 6x6 unit cells and 
/i = 0. 



from the ice state to the dodecamer ordered state is a phase transition since the translational 
symmetry is broken due to the dodecamer ordering. 

Temperature dependences of Dq and specific heat Cs for Nq = 6x6 and [i = are also 
shown in Fig. 5. Dq increases rapidly at around the sharp peak temperature in the specific 
heat. The temperature at the sharp peak is almost consistent with the inflection point of 
Dq. Thus, we conclude that the sharp peak corresponds to a transition to the dodecamer 
order driven by the effective long-range interaction due to the kinetic energy gain. In fact, 
the majority of the entropy release in the DESI system is caused by the dodecamer order as 
shown in Fig. 3(b). The similar behavior is observed in the system with 4x4 unit cells. Note 
that the sharp peak is smeared out by the finite-size effect. 
3. 3 Stability of the dodecamer ordered state 

A stability and an origin for the dodecamer ordered state have been discussed in Ref. 
11. Let us now present brief summary. (1) The dodecamer order is realized in a wide doping 
region 1/3 < n < 1/2 (—0.3 < /U < 0.2), where n is the number of particles per site. (2) The 
dodecamer order is consequence of a selection of a certain path, along which electrons can 
move easily, i.e., the kinetic energy gain. The result of (2) is concluded by the density of states: 
An energy gap is not opened by the dodecamer ordering. From this result, we concluded that 
the dodecamer ordered state is metallic. The result of (1) also means that the dodecamer order 
is not sensitive to the Fermi surface. In this way, the results of (1) and (2) are consistent with 
each other. Considering the DE mechanism which gains the kinetic energy by broadening the 
bandwidth, the weak n-dependence of the dodecamer ordered state is understood as follows. 
The ferromagnetic interaction is dominant among the localized spins at around [i = where 
the kinetic energy is mostly gained. Although there is the spin structure that is characterized 
by the finite wave vectors, the spin structure is not originated from the nesting of the Fermi 
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Fig. 6. Temperature dependence of n for the system with Nc = 6x6 and fi = 0, where black squares 
represent n at each temperature. The dashed line represents the approximate phase boundary for 
the dodecamer ordered state T* estimated from Dq [see Ref. 11 for further details]. There is the 
two-step structure which corresponds to the crossover to the ice state and the transition to the 
dodecamer ordered state as well as the entropy. 

surface, but from the frustration, i.e., the uniaxial anisotropy. 

We show the temperature dependence of n for the system with Nq = 6x6 and fj, = in 
Fig. 6. The two-step behavior can be found, which corresponds to the crossover to the ice state 
and the phase transition to the dodecamer ordered state, n decreases gradually from T ~ 1 
to 0.2. In the ice state, n keeps the almost same values which can be regarded as a plateau 
structure. Going down to further low temperatures, n shows a rapid reduction at around the 
sharp peak of C$. The behavior of n is consistent with that of the entropy. The values of 
n at [i = are also included in the region | < n < ^. It is expected that the temperature 
dependence of n is not affected drastically by changing // around zero. Thus, similar arguments 
are considered to be useful in the region of \i where the dodecamer ordered state is observed. 
3.4 Possibility of further transition in the dodecamer ordered phase 

In the dodecamer phase, e.g., Fig. 1(a), we can find that dodecamers are two- fold degener- 
ate. One has out-spins on the tips, and the other has in-spins as in Fig. 7(a). Two dodecamers 
can be regarded as "pseudo-up" and "pseudo-down", respectively. Then the dodecamer or- 
der can be described using pseudospins on a triangular lattice. Due to the DE mechanism, a 
pseudo-up and a pseudo-down are connected by the two bonds with the transfer integral t\ 
[see Fig. 7(b) (right)], while the same pseudospins are connected by the two bonds with the 
transfer integral ti [see Fig. 7(b) (left)]. There is an AF interaction between the pseudospins 
since t\ > t^. Thus, in the dodecamer phase, the problem of the DESI system on the kagome 
lattice can be mapped onto that of psuedospin systems with the effective AF interaction on 
the triangular lattice. When the random pseudospin configurations are realized, it is expected 
that the entropy S S (T) is ^log2 ~ 0.06 since the dodecamer has two degrees of freedom. 
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pseudo-up pseudo-down 
(a) two species of dodecamer and pseudospins 




(b) inter-dodecamer interaction and pseudo-AF interaction 



Fig. 7. (a) Two species of dodecamer: pseudo-up and pseudo-down. A pseudo-up and a pseudo-down 
have out-spins and in-spins on pointed tips of the dodecamer, respectively, (b) Inter-dodecamer 
interaction: An inter-dodecamer interaction behaves like the AF interaction, because there are 
two bonds with the transfer integral t\ between a pseudo-up and a pseudo-down. On the other 
hand, there are two bonds with £2 between the same pseudo-spins. Solid and dashed double lines 
represent pseudo-AF and pseudo-ferromagnetic bond, respectively. 
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Fig. 8. Temperature dependences of the specific heat Cs and the entropy S§ of localized spins for 
the system with 6x6 unit cells. The upper and the lower dashed line represent log 2 ~ 0.06 
and 0.34/12 ~ 0.03, respectively, where 0.34 is the value of the residual entropy of the AF Ising 
model on the triangular lattice. 



However, it is considered that some of the entropy Ss (T) per site are released due to the AF 
interaction. For example, in the AF Ising model on the triangular lattice, the exact residual 
entropy is known to be about 0.34 per site, 1 which corresponds to the residual entropy of 
0.34/12 ~ 0.03 per site for the present model. Magnified views of Cs and Ss at the sharp 
peak are shown in Fig. 8. From the results as in Fig. 8, it seems that 0.03 < Ss < 0.06 in the 
dodecamer phase. Thus, we may read off that there still remain degrees of freedom for local- 
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Fig. 9. The pseudo-AF phase. The uparrows and the downarrows represent the dodecamer with the 
out-spins on the tips and that with the in-spins, respectively. 

ized spins even in the dodecamer phase. From this result, disordered psuedospin states which 
satisfy the AF interaction seem to appear in the dodecamer phase. However, at further low 
temperatures, a pseudospin ordering may be expected to be realized considering the effective 
long-range nature due to the DE mechanism. A dodecamer phase without/with a long-range 
order of the pseudospins is called "pseudo-para phase" / "pseudo-AF phase", hereafter. The 
former case has only a translational symmetry breaking due to the dodecamer order, while 
the latter case also has a spin symmetry breaking due to some kind of pseudospin ordering 
in addition to the translational symmetry breaking. In Fig. 9, we show a possible pseudo-AF 
phase as a candidate for the ground state. In the pseudo-AF phase, the strong bonds with the 
transfer integral t\ linearly align in two directions, while the weak bonds with ti align in the 
other direction. The pseudo-AF phase has a rotational symmetry breaking in the lattice space 
in addition to the translational and the spin symmetry breaking. In the pseudo-AF phase, 
further kinetic energy gain may be expected because of the linear alignment of the strong 
bonds in two direction, that is, the motion of electrons on the pseudo-square lattice will be 
possible. As far as the DESI system with 4x4 and 8x8 unit cells has been investigated, 
our numerical calculations indicate that the pseudo-AF phase has the smallest energy in all 
dodecamer ordered states, i.e., the ground state of the DESI model on the kagome lattice. 
According to these considerations, we expect that there is a possibility of a phase transition in 
the dodecamer ordered state, i.e., a transition from the pseudo-para phase to the pseudo-AF 
phase. However, the ground state in the thermodynamic limit cannot be clarified so far, and 
more detailed analysis at low temperatures for larger system sizes are needed. 

4. Discussion and Summary 

In the DESI system on the kagome lattice, the long-range interaction due to the DE 
mechanism plays an important role to stabilize the dodecamer ordered state. In fact, a similar 
situation has been seen in the spin ice system. 13 The spin ice system has been discussed by two 
models: One has a ferromagnetic interaction due to the n.n. Ising interaction, and the other 
has that due to the long-range dipolar interaction. In particular, the latter model is called a 
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(a) hexamer 

Fig. 10. (a) Hexamer in ZnCr204. The S = 3/2 localized spins align antiferromagnetically in an 
arbitrary direction, (b) Octamers in CULY2O4: (left) Octamer composed of eight Ir 3+ (S = 0) 
which is represented by open circles, (right) Octamer composed of eight Ir 4+ (S = 1/2) which is 
represented by open circles, where the double lines represent the dimer-singlets of Ir 4+ . 

dipolar spin ice model. According to MC calculations in the dipolar spin ice model, a double- 
peak structure in a specific heat has been observed as in the DESI system: A broad peak at 
higher temperature and a sharp peak at lower temperature. 15 The sharp peak corresponds 
to the spin cluster ordered state. Such a cluster ordered state is not observed in the spin 
ice model with the n.n. Ising interaction. These results indicate the long-range nature of the 
interaction plays an important role to create a cluster ordering even in the frustrated spin 
systems as well as the electron systems. Although the spin ice system is an insulator, the 
DESI system is considered to be metallic. In this way, a cluster ordering might be generic to 
frustrated systems. In general, the frustration suppresses the usual ordered states and leads to 
the degenerate states at low temperatures. In the frustrated systems, the residual interaction 
selects the cluster ordered state among the macroscopic degenerate states. 

Let us consider a possibility of the cluster order in realistic systems with the frustration. 
As mentioned in the introduction, there are several examples of a cluster order in nature: 
the hexamer in ZnC^CU 9 [see Fig. 10(a)], the octamer in CUIX2S4 10 [see Fig. 10(b)] and so 
forth. In the cubic spinel ZnCr204, Cr 3+ (S = 3/2) resides on the pyrochlore structure and 
interact with each other antiferromagnetically. 16 At low temperatures, a hexagonal AF spin 
cluster formation is realized and the pyrochlore sites are tiled by the hexamers. 9 In a thiospinel 
CuIr 2 S4, mixed-valence states of Ir 3+ (S = 0) and Ir 4+ (S = 1/2) are realized in the insulating 
phase at low temperatures and they occupy octahedral sites. 17-19 In this case, the system have 
a charge frustration as well as a spin frustration. In the insulating phase, Ir sites are tiled by 
two kinds of octamers: One consists of eight Ir 3+ and the other consists of eight Ir 4+ with four 
dimer-singlets. 10 In this way, the cluster order can be realized in realistic frustrated systems. 
Such a cluster order may be explained by including the effect of the residual interactions which 
are neglected in the non-frustrated system. 

In conclusion, there are at least two steps for the entropy release in the DESI system on the 
kagome lattice, which is indicated by the peak structures in the specific heat. The ice state 




(b) octamer 



13/18 



J. Phys. Soc. Jpn. Full Paper 

appears at intermediate temperatures due to the n.n. ferromagnetic interaction. From the 
estimated entropy, it has been clarified that there is a large number of thermally degenerate 
states in the ice state. As the temperature is lowered, a large part of the entropy in the 
system is released by the dodecamer ordering with the translational symmetry breaking. The 
dodecamer order is driven by the kinetic energy gain and the frustration: The dodecamer is 
the kinetic-energetically stable path for electrons. Furthermore, a possibility of a transition in 
the dodecamer phase is indicated by our MC results: (1) The estimated entropy Ss seems to 
remain in the dodecamer ordered phase because of the existence of two kinds of dodecamers 
which can be regarded as the pseudospins. (2) The pseudo-AF phase is the ground state in 
the system with small-size clusters. These result possibly indicate that there is a transition 
from the pseudo-para phase to the pseudo-AF phase to gain the further kinetic energy. It is 
expected that the entropy in the system will be completely released by the ordering of the 
pseudospins, i.e., two kinds of dodecamers. However, more detailed analysis are needed to 
clarify this point. 
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Appendix: Monte Carlo Method 

Here, we explain the application of the MC method for the Hamiltonian (5) in detail. 20 
We define abbreviations {cj,c|} = (c±,--- ,cn s ,c\, ■ ■ ■ , c N ) and {t.i} = (ri,--- ,tn s ) which 
represent degrees of freedom for electrons and localized spins (adiabatic fields), respectively. 
Using Eq. (5), the partition function of the system in the grand canonical ensemble is written 
as 



Z = TrcTrp exp 



TrcTrp exp 



-P\Hk 



* (n, Tj) \ c\cj + /i.e.) - n 

i 



(A-l) 



(A-2) 



where Trp and Trc represent traces over {cj,c|} and {t^}, respectively. And, \i represents 
the chemical potential. Here, we define the kinetic energy part of the Hamiltonian as = 
j) * ( Ti ' T ^ ( c l c i + ^- c -)- O nce a localized spin configuration {r^} is given, it is easy to 
diagonalize the Hamiltonian (5) because it is represented by the quadratic form of fermion 
creation and annihilation operators, i.e., As x As matrix: Hk = c l^ij c j^ where Hij = 

t(ri, Tj) for n.n. pairs of i and j, otherwise 0. Thus, the trace operation Trp is easily performed 
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through the energy eigenvalues £fc({Tj}) of for a given configuration {r^}. Finally, the 
partition function (A-2) can be written as 

Z = Tr c exp[-/3F cff ({r i })], (A-3) 

with the effective free energy 



Fesd-ri}) = ~k B T log (1 + exp {-P [e k ({Ti}) - »]}) . 
k=l 

Probability distribution of a certain configuration {-n;} is given by 

P({n})=exp [-/3F eff ({ri})] /Z, 



(A-4) 



(A-5) 



whose numerator is regarded as the Boltzmann factor. 

Thermal average of arbitrary observable A({cj, c}}, {t^}) is given as follows in the grand 
canonical ensemble: 



Tr c Tr F 


~A({ Cl ,4},{n})e-^~ 


Tr c Tr F 


\-f3H~ 





(A-6) 



Using Eq. (A-5), the thermal average is also written as 

(M{ci,4h{Ti})) = ((i({Q,4},{r i })) F ) C; (A-7) 
where (• • • )f represents the grand canonical average for a fixed configuration {rj} defined by 



<■■■>! 



Tr F 




Tr F 







(A* 



while (■ ■ ■ )c is the thermal average over configurations of {t^} defined by 

(•..)c = Trc[-P({Ti})]- (A-9) 

Note that, in the thermal average for observables without the electron operator {cj,c|}, e.g., 
vf (Eq. (7)) and di (Eq. (10)), (• • • ) becomes (• • • )c- 

The specific heat is defined as a fluctuation of the energy in the form 

(H 2 ) - (H) 2 



C(T) 



rp2 



(A10) 



C(T) can be decomposed into electronic part Ce(T) and spin part Cs(T) as follows: 

C(T) = C E (T)+C S (T), (A-ll) 

where 

(A-12) 



C E (T) = ^(E 2 -E 2 ) C , 
C S (T) = -L ((E 2 )c - <£i) 2 c ) • 



(A-13) 



15/18 



J. Phys. Soc. Jpn. 



Full Paper 




CO 



Fig. A-l. Temperature dependence of Ce and Se- Ce has a peak at T ~ 0.6, which is almost the 
same order of t\ and £2- 

Here, we set E\ = (H)f and E2 = (H 2 )f. Eq. (A- 12) is given by the electronic energy fluctu- 
ation averaged over configurations {t^}, while Eq. (A-13) is given by the energy fluctuations 
among spin configurations {Tj}. Thus, it is easily found that Eq. (A- 12) does not include fluc- 
tuation of localized spins while Eq. (A-13) does not include fluctuation of electrons. E\ and 
E2 are calculated as 

E 1 =Y,^({n})f F (e k ({n})), (A-14) 

k 

E 2 = Y / e k ({n}) 2 Me k ({r l })) 

k 

+ X)e*({ri})ei({T i })/i.(e fc ({r i }))ys.(e,({r i })). (A-15) 

k^l 

where fp(x) is the Fermi distribution function given in the form, 

In the grand canonical ensemble, the entropy of electronic part per site in the DESI system, 
Se(T = 00) is log 2 as follows: 

S E (T = 00) = 

= (log2) c 

= log2. (A- 17) 

In the DESI system, the entropy of localized spin is log 2, since each localized spin has two 
degrees of freedom. Therefore, total entropy per site in the DESI system is log 4. Temperature 
dependences of Ce and Se for the system with 6x6 unit cells are shown in Fig. A-l. A peak 
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can be found at around T = t\ ~ £2 in Fig- A-l. In the electronic part of the specific heat, Ce, 
any remarkable structures are not observed at sufficiently low temperatures T ~ 0.01 where 
the dodecamer order is realized. 
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